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ABSTRACT 


An experimental apparatus which will allow the study of 
the influence of ground effects on a heated cylinder placed 
in a crossflow of air was designed and constructed. 

The effects of a wall on the heat transfer character- 
istics of a heated cylinder were investigated at a Reynolds 
number of 153,000 for cylinder to wall spacings ranging from 
d/r = 0.0 to d/r = 5.33 (maximum), where d = gap width, and 
r = cylinder radius. Pressure coefficients were obtained at 
five degree intervals around the cylinder for each experiment. 

Experimental results obtained for d/r = 0.0 agreed well 
with previously reported investigations of a cylinder 
attached to a wall. Results obtained for d/r = 5.33 agreed 
well with previous work done with free cylinders. Experiments 
conducted at intermediate values of d/r resulted in minimum 
and maximum average Nusselt numbers, suggesting that gap 
Spacing may exert a strong influence on the behavior of the 


Flow field. 
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Pee ROUUCTION 


The principal objective of this work has been to design, 
construct and test an experimental apparatus which will allow 
the study of the influence of ground effects on a heated 
cylinder placed in a crossflow of air. Previous work of this 
nature has concentrated primarily on free cylinders. 
Recently, McComas (1) investigated the pressure and heat 
transfer characteristics of a heated cylinder attached to a 
wall. McComas determined the local and average Nusselt 
number, and the pressure coefficient, on the surface of the 
cylinder for Reynolds numbers varying from 53,000 to 153,000. 
This work is a follow-on study of the project initiated by 
McComas, differing in that the apparatus has been designed 
and constructed such that the plane surface can be varied 
from the attached position to approximately 3 cylinder 
diameters away. Due to the limitations of the wind tunne! 
used by McComas, (described in detail in (1)), hotwire 
anemometry and flow visualization techniques were not able 
to be used in his investigations. The present experimental 
apparatus is located in the Mechanical Engineering Department 
wind tunnel in Building 500 on the Naval Postgraduate School 
campus. This wind tunnel, described in detail in section III, 
has a speed range of 10 to 175 miles per hour, and is equipped 
with a traversing mechanism and hotwire anemometer which a 


be used to obtain velocity information throughout the test 





section. The tunnel also has a capability for flow visual- 
ization using smoke filaments. 

| Experiments were conducted at a free stream Reynolds 
mmper of 153,000 for cylinder toplate spacings of d/r = 0.0, 
oze>, 0.5, 1.0, 2.0, 4.0, 5.33 (maximum separation) where 

d = gap width and r = culinder radius (Figure 1). The flow 
field in the vicinity of the cylinder was probed using the 
hot wire anemometer. The velocity distribution suggested 
that Flow blockage was a factor in the experiment, and that 
the Reynolds number should be based on a corrected free 
stream velocity. The correction produced an effective 
Reynolds number of 161,006. Using this corrected Reynolds 
number, the data taken at d/r = 0.0 agrees with the results 
of McComas, and the data at d/r = 5.33 agrees with the work 
oamraeld (2). 

As was the case with McComas, experiments with the 
cylinder attached to the plateshowed a critical flow pattern 
to exist over the top half of the cylinder. As the cylinder 
to plate spacing was increased, this flow pattern changed 
imegecyitical to subcritical over a narrow range of plate 
spacing (at about d/r = 0.09). The flow over the bottom 
half of the cylinder exhibited critical flow as soon as a 
Space existed. This flow also changed from critical to 
subcritical, but the change was less pronounced and it was at 


a larger cylinder to plate spacing (at about d/r = 1.2). 
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Figure 1. Schematic diagram showing determination of d/r. 
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Another phenomenon which was seen to occur with increasing 
d/r was the gradual shift of the stagnation and laminar 
Separation points over a range of 15 to 20 degrees. 

In calculating and plotting the average Nusselt number 
as a function of gap spacing, an interesting and unexpected 
trend was revealed. It was found that the minimum average 
Nusselt number did not occur when the plate was attached 
momune cylinder, but at a d/r = 0.25. Further, this curve 
also indicated a maximum point at d/r = 2.0. Although the 
current data are incomplete, (further investigation is 
needed over a range of Reynolds numbers) the Nu versus d/r 
curve at Re = 153,000 suggests that gap spacing may exert a 
strong influence on the behavior of the flow field. Future 
studies employing both flow visualization as well as heat 
transfer measurements would aid in establishing a plausible 
explanation for the trends observed in the curve of average 


Nusselt numbers. 





ee DA OChGROUiaD 


Although a great deal of investigation has been devoted 
to the flow and heat transfer characteristics of cylinders 
mar from boundaries, see Meyer (5) and Field (2) for a list 
of references, virtually no references have been found 
relating to the heat transfer and flow characteristics of a 
cylinder placed in close proximity to a plane surface. The 
ferew ical applications of such information would extend from 
aiding in conducting thermal analysis of container stored 
ordnance, rockets and missiles to the determination of the 
optimum size and spacing for heat exchanger tubes. 

Although McComas (1) experimentally determined local 
and average heat transfer coefficients on the surface of a 
uniformly heated cylinder attached to a wall, McComas found 
no references to eitner theoretical or experimental work on 
this subject. McComas proposed a theory, based on theories 
that had been postulated for forward and rear facing steps, 
that quantitatively explained the trends in his heat transfer 
data. Other than McComas’ work, no further references have 
been found dealing with attachedcylinders. 

Since the fluid flow field and thermal field are 
intimately coupled in a forced convection process, an under- 
Standing of the flow fieid is mandatory if one is to under- 
Stand the heat transfer process. A thorough discussion of 


miemeriia flow and heat transfer characteristics of a free 





cylinder can be found in (3), and a proposed but as yet 
unverified description of the flow around an attached 


cylinder can be found in (1). 
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Pies tCeR Nien ALeAreARAIUS 


A. WIND TUNNEL 

All tests were conducted in a subsonic, open circuit, 
wind tunnel manufactured to Navy specifications by the 
Kenny Engineering Corp., of Monrovia, California (Figure 2). 
The overall length of the tunnel was 40 feet, with an 
entrance bell 5 feet, 5 inches wide by 7 feet, 5 inches high. 
The test section was 20 inches wide by 28 inches high at 
the entrance, expanding to 3] inches high to allow for 
boundary layer growth. These dimensions resulted in a 
contraction ratio, based on area, of 10:1. The tunnel was 
powered by a 75 HP, 1750 rpm variable speed motor driving a 
45 inch diameter adjustable pitch vaneaxial fan capable of 
supplying air at 10,000 to 70,000 cubic feet per minute. The 
tunnel was equipped with six graded mesh screens to control 
turbulence and a 24 inch manometer tube tapped into the 
stilling chamber and test section for velocity measurement. 
The tunnel was designed to operate at velocities ranging 
eon )0 moh (14.67 ft/sec) to 175 mph (256.67 ft/sec). 
Turbulence intensity was taken using a hot wire anemometer 
and found to be approximately 0.2% at a free steam velocity 


of 98 ft/sec (see Appendix D). 
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B. HEAT TRANSFER EXPERIMENT 

The test cylinder was essentially three cylinders 
attached end to end (see Figure 3). The two end sections 
were hollow softwood cylinders supported by a hardwood 
brace, curved to match the cylinder contour. The central 
cylinder was a Sheet of resistive paper (Temsheet) fitted 
into stepped portions of the end cylinders, and attached to 
the end cylinders and hardwood brace by double sided adhesive 
tape. The three cylinders then formed one cylinder 20 inches 
long, to correspond to the test section width, and having a 
diameter of 3 inches, chosen to allow maximum resolution of 
the liquid crystals without excessive blockage within the 
test section. The assembled test cylinder was tightly 
packed with glass wool to minimize heat transfer into the 
cylinder, and thereby reduce the time needed to reach steady 
state after a voltage was applied to the paper. 

The resistive paper used was a commercially available, 
carbon impregnated paper, Temsheet, with a nominal thickness 
om 0.039 renee: The width was sized to conform to the 
Stepped circumference of the end cylinders (9.422 inches) and 
the length was 12 inches. The circumference of the resistive 
paper cylinder was marked in 5 degree increments, in two 
locations 2 inches apart. This was the test section used in 


all heat transfer experiments. The entire outer surface of 
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the resistive paper was covered with three coats of NCR S-43 

liquid crystals having event temperatures of: 

red: 109. 8°F 

green: 111.0°F 
mi2z.3 F 


blue: (Field (2)) 


Power was applied to the paper by means of wires lead 
inside the cylinder through one end. The wires were attached 
to aluminum foil tape by means of alligator clips, and the 
edges of the tape were painted with silver paint to assure 
good electrical contact with the paper. The nominal distance 
between the tape/paint was 4.8125 inches, resulting in a 
surface area for power dissipation of 45.34 in’. 

The power source used was a Lambda regulated power 
supply, model LK 345A FM. The voltage applied by the Lambda 
was measured by a Fluke 8100A digital multimeter. This 
multimeter was also used to measure the resistance of the 
resistive paper. Air stream temperature was sensed by a 
copper-constantan thermocouple installed in the stilling 
chamber. The thermocouple output was read on a Leeds & 
Northrup millivolt potentiometer. 

The "wall" used in these experiments was manufactured 
from 1/2 inch plexiglass. This platewas 46 inches long, the 
maXimum practica! length, to minimize the influence of the 
trailing edge on the flow patterns around the cylinder. The 
width of the plate was 20 inches, to conform to the test 
section width, and it was beveled on the under side of the 


leading edge at an angle of approximately 12 degrees. 


ia 





The plate and cylinder were placed in the tunnel 
horizontally, and in such a way that the leading edge of the 
plate extended four cylinder diameters (12 inches) ahead of 
the cylinder centerline. The cylinder was attached to the 
test section walls by portions of the end cylinders pro- 
truding through holes cut for that purpose. The plate was 
capable of being moved vertically by means of a threaded 
rod extending through the tunnel floor, and a rear brace 
which locked in position by means of a split nut (Figure 4). 
The plate was capable of being located at any position 
relative to the cylinder, from attached, to a maximum 


possible distance of 8 inches away (Figure 1). 


on Bae URE CURT LOLEN | EXPERIMENT 

The experiments conducted to find the variation of 
pressure around the circumference of the cylinder used the 
previously described plexiglass plate as the wall and a 
3 inch diameter aluminum cylinder fitted with four pressure 
taps 90 degrees apart (Figure 5). The circumference of the 
cylinder was scribed in 5 degree increments. Plastic tubing 
led from the taps through one end of the cylinder to a 
manometer bank. 

The manometer bank consisted of several tubes that were 
inclined at an angle of 30 degrees from the vertical. All 
four cylinder pressure taps as well as the stilling chamber 
and test section pressures were read from this manometer 


bank. 
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Photograph of the cylinder used to collect pressure data. 
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[i oeeeAPERIMENTAL PROCEDURE 


The present work was conducted utilizing different test 
facilities than those used in previous investigations. The 
wind tunnel was smaller, and of the open circuit variety. 
Due to experimental considerations, the cylinder size and 
iemctation were different, and the problem of varying the 
relative separation of the plate and cylinder was unique to 
this work. In order to assure that meaningful and accurate 
data would be obtained, it was necessary to compare initial 
results with previous work. Since both Field and McComas 
performed experiments at a Reynolds Number of 153,000, it was 
decided to check the present apparatus at that Reynolds 
Number. For this reason, an air stream speed of 98 ft/sec 
was used for comparison experiments. Because of time 
limitations, all of the data for the present work were taken 
eumeo 1ce/sSec, corresponding to a Reynolds Number of 153,000, 


where Reynolds number is defined as: 





Ud 
Re = 
Vv 
Re = Reynolds number 
i icccmormediimyelOCltyede testesection 
entrance, ft/sec 
d = cylinder diameter, ft 
v = kinematic viscosity of air, ft°/sec 
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A. HEAT TRANSFER EXPERIMENT 

The same experimental procedure used by Field and 
“McComas was used in the present experiments. For complete- 
ness, a description of the procedures will be repeated here. 

Prior to the start of the experiment a sufficient 
voltage was applied to the cylinder to cause the liquid 
crystals to reach their event temperature and change color. 
This voltage was usually from 12 to 15 volts. After about 
a half hour, time to allow the cylinder and glass wool 
packing to come to a steady state temperature, the power 
supply was turned off and the multimeter was switched from 
DC voltmeter to ohmmeter. In this manner an accurate value 
of the paper resistance was determined at a temperature which 
corresponded to the liquid crystal event temperature. The 
multimeter was then switched back to DC voltmeter and the 
wind tunnel was started and brought up to speed. 

ies wind tunnel velocity range and test section velocity 
profile were examined prior to obtaining experimental heat 
transfer data (see Appendix D). It was determined that the 
velocity profile was uniform at the test section entrance, 
and that the turbulence intensity was approximately 0.22. 
After numerous checks with the hot wire anemometer, it was 
decided to use the installed wall manometer to set the tunnel 


speed of 98 ft/sec. 
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After the tunnel reached operating speed, the voltage 
was adjusted such that the liquid crystal color bands would 
Prpear at various angular positions around the cylinder. 

The first color band to appear would be at the lowest local 
Nusselt number on the cylinder. For tests with a d/r greater 
then 0.25 (where d = distance from plate to cylinder rim, and 
r = cylinder radius) this was always the laminar separation 
point (Figures 13-19). As voltage was increased, the color 
lines would appear in one, or as many as 5 or 6, angular 
locations all having equal Nusselt numbers. In this manner 
the cylinder surface temperature could be compared to the 

air stream temperature as determined by a thermocouple 
mounted in the wind tunnel stilling chamber. Using this 
temperature difference, and calculating the power dissipated 
from the voltage applied, the local heat transfer coefficient 
could be calculated for the circumference of the cylinder. 
Because of the way the cylinder was constructed, readings 
within 15 degrees of the seam were considered unreliable. 

The majority of Nusselt number data were collected with the 
seam located at 180 degrees. The Nusselt number in the 
vicinity of 180 degrees was determined by turning the cylinder 
such that the seam was facing forward at 0 degrees. Several 
experiments were also made with the seam at positions of 

135, 225, and 315 degrees. The final position of 0 degrees 
was chosen for the seam after it was observed tnat the value 


and angular location of the Nusselt number at the laminar 
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Separation point, and in the vicinity of 130 and 210 degrees, 
was consistent with the values obtained with the seam at 
| 180 degrees. 

For all experiments, distinct, repeatable values of the 
Nusselt number versus angular location were obtained. The 
gradient as determined by the color bands, was sharp, with 
the exception of a region near the forward stagnation point, 
which appeared nearly isothermal. The gradient was also 
indistinct in the turbulent boundary layer, where the 
"scrubbing" action of cooler air caused the Nusselt number 
to change with time. Therefore, the Nusselt number shown 
for the turbulent regions must be considered to be an average 


value. 


B. Peeosvinem COEFFICIENT EXPERIMENT 

With the pressure cylinder mounted in the wind tunnel, 
the manometer bank was leveled, pressure leads checked for 
tightness, and the wind tunnel brought up to a speed of 
98 ft/sec, corresponding to a Reynolds number of 153,000. 
The cylinder was positioned such that the pressure taps were 
initially at 0, 90, 180, and 270 degrees. The manometer 
readings for the four taps, the stilling chamber, and the 
test section were recorded. The cylinder was rotated in 
5 degree increments until pressure readings were obtained for 


the entire circumference. 
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The pressure coefficient was determined as follows: 


where: 








pressure coefficient 

density of air 

air stream velocity (test section entrance) 
pressure on cylinder surface 

stagnation pressure (stilling chamber) 


air stream static pressure (test section 
entrance) 


ed 





Vee Reo LS 


Initial heat transfer and pressure coefficient experi- 
ments were made with the plate and cylinder attached 
Mayr = 0.0)" A thin felt strip was placed between the two 
to assure a tight seal. The data obtained in this configura- 
tion agreed well with the results of McComas (see Figure 6). 
It can be seen that the trends in the data are identical, 
that local maximum and minimum values are within experimental 
uncertainty (see Appendix B for uncertainty analysis), and 
that a local maximum and local minimum were again observed 
near the point of attachment of the plate and cylinder. This 
is further evidence of the trapped vortex phenomenon as 
proposed by McComas. 

The next series of tests were made with the plate as 
far removed from the cylinder as possible, d/r = 5.33, or 
d = 8 inches. As there appeared to be some blockage, the 
method of Pope (4), as used by Meyer (5), was used to 


determine the necessary correction: 


R = (ae + e)Re_ 


“corr 
where: 

eee = be “wb 
and oan is the solid blockage correction factor: 


ee ne cry? 
Sb or it 
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and eat is the wake blockage correction factor: 


eee l/c (D/H) Ch 


where: 
D = cylinder diameter, inches 
yr = cylinder radius, inches 
H = total tunnel height, inches 
Ch Mialg “Coerficient (see Meyer p. 107) 
Therefore: 
2 55.2 1 ee : 
and: 
aaa = e055) 1535000 
No ain = 161,000 
Using the on to compare the results of this experi- 


ment with the work of Field, it was found that the results 
compare favorably with Field and with the theory of Schuh (6) 
(see Figure 7). It should be noted that Field experienced 
critical flow at Re = 149,000, while the experiments made in 
this work were exclusively subcritical for d/r = 5.33. This 
phenomenon could be attributed to higher free stream turbulence 
or cylinder deformation, either of which could have induced 


critical flow in Field's experiments. 
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Figure 7. Comparison of the heat transfer results 
fOmcd/s~we= 5.2.33 and See = 161,000 with 


tnemresurtts of Field at Re = 148,000 
and with the theory of Schuh 
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Figures 13 through 19 are polar plots of the local 
Nusselt number and pressure coefficient for values of d/r 
ranging from 0.0 to 5.33. In examining the seven plots for 
Significant trends or unusual phenomena, two items of 
interest were immediately apparent: one was the change from 
critical to subcritical flow, the other was a shift in the 
location of the stagnation point and the laminar separation 
points as the plate was moved away from the cylinder 
(see Figure 8). 

ice tLnensitNOnetronecriti1cal to subcritical flow was 
investigated in order to ascertain the d/r ratio associated 
with the transition. As the flow was critical with the 
plate attached, and subcritical over the top at d/r = 0.25, 
the first transition point of interest was located in the 
mmeerval 0.0 < d/r < 0.25. By alternately moving the plate 
and adjusting the voltage, the transition point from critical 
to subcritical flow over the top of the cylinder was found 
to be approximately 0.08 < d/r < 0.10. 

As the flow over the bottom of the cylinder underwent 
a transition from critical to subcritical flow in the interval 
1.0 < d/r < 1.5, this region was investigated further. Using 
the same procedure of carefully adjusting spacing and voltage, 
the transition was found to occur at approximately 


1.17 < d/r < 1.25. 
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It should be noted that while both flows described 
above are definitely critical, there is significant 
difference in the degree of criticality on the top and 
bottom of the cylinder. This can be seen from the magnitudes 
of the Nu on the data shown in Figures 9 and 10. 

In order to glean further information from the data, 
the average Nusselt number, Nu, was computed from the 


equation: 


NO 


Om 
——~ | 
u = | Nu dé 
T 
0 


which was numerically integrated by computer using Simpson's 
Rule of Quadrature with a step size of 5 degrees. The 
results are shown in Figure lla. 

It is interesting to note the unexpected shape of the 
curve. Nu drops to a minimum value, then increases to a 
maximum, and then gradually decreases to a nearly constant 
value as d/r approaches relative infinity. These results 
prompted further investigation. 

In order to help determine the influence of the local 
Nu on Nu a plot was made of Nu, for the forward (-90 to 
90 degrees) and rear (90 to -90 degrees) portions of the 
cylinder. The results are shown in Figure 11b and exhibit 
the same general shape as Nu. The drop in Nu from d/r = 0.0 
to 0.25 is explained from this plot when the sharp drop in 


ee imei yonusot tie cylinder from d/r = 0.0 to 0.25 is 
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Figure 11b. Maximum Nusselt numbers on the front and 
and back of the cylinder versus d/r. 
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noted. (The maximum local Nu for d/r = 0.0 occurs in the 
separation bubble, in the vicinity of 90 degrees, and is not 
) 


Further investigation of the Nu curve lead to an attempt 


used in the plot of NIU eas 
to see how the voltage applied to the cylinder would vary if 
one location on the cylinder was maintained at a constant 
temperature as the plate spacing was varied. The locations 
150 to 165 degrees were chosen for this experiment. For each 
d/r the voltage was adjusted until a distinct color band 
appeared at these locations. Plotting the voltage applied 
versus d/r resulted in a curve of similar shape to Nu and 
Ue a x oom eiaure: 2). | 

In order to verify the curvature of the Nu plot, a 
complete set of data was taken at d/r = 3.0. This data is 
not presented on a separate polar plot but Nu was calculated 
and the value obtained enabled the shape of the Nu curve to 
be verified. 

A possible explanation for the observed variation in 
Nu could be related to the amount of "cold" free stream air 
which has access to the cylinder surface. At d/r = 0.0, the 
flow is critical due to the fact that all of the air approach- 
,ing the cylinder is forced to go "up and over." As soon as 
the plate is separated the flow has two paths. As soon as 
Space is sufficient to provide a favorable path without large 
pressure gradients, enough of the air approaching the cylinder 


passes beneath to cause the flow over the top to transition 


So 
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to subcritical. Now the air passing below the cylinder is 
experiencing critical flow, possibly due to the nozzle 

effect produced by the rapidly decreasing flow area below the 
cylinder. The heat transfer in this region does not increase 
due to the fact that much air has already been heated passing 
over the front of the cylinder (i.e., the thermal and 
hydrodynamic boundary layers are a significant fraction of 
the gap size). 

As d/r is increased beyond 0.25, the larger opening 
allows more cold air, outside the boundary layer, to approach 
the bottom of the cylinder. This cold air carries off some 
heat by conduction and also mixes with the warm air in the 
separation bubble and turbulent boundary layer, resulting in 
increasing Nu. 

The gradual decrease in Nu froma high value at 
d/r = 2.0 is not understood at this time. It was expected 
that Nu would approach some constant value as the influence 
of the plate was reduced. It can be seen from Figure lla 
that the uncertainty band nearly encompasses the decrease 


in Nu. Although the behavior of the Nu, and voltage curves 


ax 
Suggest that the trend may be genuine, further investigation 
-is needed to determine if Nu actually has a peak, or in fact 
approaches a constant value from some minimum point at a 


eiec |! d/r. 
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Figure 13. Polar plot of Nusselt number and pressure 
coefficient versus angle for d/r = 0.0. 
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Figure 14. Polar plot of Nusselt number and pressure 
coefficient versus angle for d/r = 0.25. 
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meoure '5. Polar plot of Nusselt number and pressure 
coefficient versus angle for d/r = 0.5. 
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Polar plot of Nusselt number and pressure 
coefficient versus angle for d/r = 1.0. 
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Figure 17. Polar plot of Nusselt number and pressure 
coefficient versus angle for d/r = 2.0. 


45 





Figure 18. 
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Polar plot of Nusselt number and pressure 
coefficient versus angle for d/r = 4.0. 
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muoure 19°" Polar plot of Nusselt number and pressure 
coefficient versus angle for d/r = 5.33. 


47 





a 


VI. CONCLUSIONS AND RECOMMENDATIONS 


ine phitdrysonujective of this work, to design, construct 
and test an experimental apparatus which will allow the study 
of ground effects on a heated cylinder placed in a cross flow 
of air, has been achieved. The comparisons made with 
previous work are within experimental uncertainty, and 
indicate that the results of further experiments made with 
this apparatus can be relied upon. The liquid crystal 
technique for heat transfer analysis has proved to be an 
excellent method for determining quantitative heat transfer 
data. The results of the one plate geometry and Reynolds 
number investigated in this work are shown in the plots of 
Nu and C, versus theta. Figure 8 shows the smooth shift 
in the location of key points which characterize flow 
phenomena (stagnation point, and laminar separation points) 
and makes it possible to approximate the value of Nu for 
configurations other than those investigated. 

Future work in this area should include a further 
investigation of the variation of Nu with increasing values 
of d/r. The first step in this investigation could be 
further experimentation with the configuration described in 
this report at various Reynolds numbers. Additional work 
Should also investigate the effect of the plate length on 
flow patterns. As the wind tunnel used in this work has a 
capability for flow visualization, a series of experiments to 


determine flow patterns visually would prove to be valuable. 
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APPENDIX A 
HEAT TRANSFER AND PRESSURE DATA REDUCTION 


The equations and procedures used to reduce the heat 
transfer and pressure data are shown below. The results of 
each experiment are presented on polar plots of Nu and Cp 


versus Theta in section \V. 


A. HEAT TRANSFER EXPERIMENT 

The relationship between power dissipated (by convection 
and radiation heat transfer), and power input by Joulean 
heating is the same as developed by Field (2) and used by 
Field and McComas (1). 


ye 
h. = 
t RA(T-T_) (0.293) 





where: hy = total heat transfer coefficient ~~ y—— 
biotite. @- 
V = voltage applied,volts 
R = total resistance of Temsheet, Ohms 
A = total heated surface area of Tewslhaats k= 
T = surface temperature of the cylinder as 
determined by liquid crystals, °F 
(ee escmmeamecemperatlife,,  F 
0.293 watt-hr _ 746 watts 


BT ~ 2545 BTU/hr 
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The Nusselt number is defined as: 
h .d 
oe Ke 
where: 


Nu = Nusselt number 
Bru 


ie = convection heat transfer coefficient, ——~—y— 
ie fi Oe 
d = cylinder diameter, ft 
7 BTU : a 
kee 0 O15 hy ft°r (air at 75 ) 


Since Nu is defined by the convective heat transfer 
coefficient, it is necessary to subtract the radiation effects 


from the total heat transfer coefficient 


This correction is made using the equivalent radiation heat 


transfer coefficient has where ee is defined as: 


2, 2 2 
ecmacde: 1.) (1 + 1.) 
where: 
h = equivalent radiation heat transfer 
" BTU 
coefficient 5 
fier tO F 
Oo = 0.1714 x 107° x A (Stefan-Boltzman Const) 
, MP ie 
C-— 


oH, > UT mpeaiacloneexchange factor between the 
cylinder and the surroundings 


T = cylinder surface temperature; °R 


T = air stream temperature, °R 
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A sample calculation for data obtained at a Re = 153,000 


and d/r = 


Nu = 


Nu = 


0.50 is shown below: 


acc we etoo VOCS 1 = 111.0°F, I = 63.5°F 


Z 2 


24.65 volt 
(13.2n) (0.314: ft°)(111-63.5°F)(.293Watt_nry 


BTU 


BTU 


10.5 a 
Mined tale 


8 el 


(Opelika 10. 5 z)(.9)(571 eee kt) 
ieee R 


e 4 523.5°R°) 


(o7] 


BTU 


hr Fto°F 


(9.5 —2,—)(.25 ft) 


fe gee ees 


pa 
BTU 


(0180 FE SEF 
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D. eee COME ICTENT EXPERIMENT 


As described previously, the pressure coefficient is 


defined as: 


i 





C. = pressure coefficient 

> = density of air, 1b/ft> 

as air stream velocity, ft/sec 

P = pressure on cylinder surface, in, 0 
2 

an = stagnation pressure, 'MH0 

p 


= air stream pressure, in 
H,0 


A sample calculation for data obtained at a Re = 153,000 


and d/r = 0.50 is shown below: 


eae Oem Pe = 6.0 P_. = 1.86 


C = - 0.99 


Se 





APPENDIX B 
UNCERTAINTY ANALYSIS 


Because of the uncertainties in the measured variables, 
jt was necessary to perform an uncertainty analysis in 
order to ascertain the degree of uncertainty in the experi- 
mental results. The method of Kline and McClintock (7) 
was used. 


The following estimated uncertainties existed: 


voltage, V + 0.1 volts C2051) 
resistance, R + 0.2 ohms (20:1) 
temperature AT + 1.0°F (20 cll ) 
area Ave 0.0 | rt° (e< 1} 
diameter dmaeo-00l ft (20:1) 
radiation 
exchange factor Ot + 0305 20:71) 


+ 


wall manometer H OOo | C2021) 


manometer bank oO Sean Pe 1) 


Total Heat Transfer Coefficient, hy: 


As shown in Appendix A, the equation used to calculate 


the total heat transfer coefficient was: 


(since the conversion from watts to BTUs is a constant, 


it does not enter into the uncertainty analysis). 


a3 





t 


For data obtained at a Reynolds number of 153,000 


end ad/r = 0.50: 


V = 24.68 + 0.1 volts 
R = 13.2 + 0.3 ohms 


AT = 47.5 + 1.0 °F 
2 


A = 0.315 + 0.001 ft 
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In some regions where T = + 2.0 °F, then: 


U) 
he. " 
Dy OOS s -or Wp, Oy s 


h = 10.5 + 0.5 


All measured quantities were read with the uncertainties 
shown above, with the result that the calculated value of 
total heat transfer coefficient had a relatively constant 


value around the circumference of the cylinder. 
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Since the angular location of the color bands does not 
enter into the calculation of Hi» the angular uncertainty 
has not been included here. This uncertainty is estimated 
Meme + 5 Gegrees,. 


Equivalent Radiation Heat Transfer Coefficient, 1 


The equation used to calculate He was: 
h, = o@(T + T_)(T* + TS) 


deeining 1, = 1/2 (T + T,) 


OW 


h = aca T 


“ 
= 


since o 1S a constant, it does not enter into the uncertain- 
ty analysis. 
For the same data used to perform the uncertainty 


analysis for h the uncertainty in h. was found to be: 


| eas 
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Nusselt Number, Nu: 


The equation used to calculate Nu was: 





Since k is a property of air, and its uncertainty is 


negligible: 





Nu = 158 + 5 


U) 
While this analysis has shown that _ SOS. and 


meet Ss degrees, it should be noted that this may be 
misleading. Because there are portions of the cylinder which 
exhibit different flow phenomena, there are different types 
of color bands exhibited. The isothermal region has no 
distinct bands, and the mottled colors which appear cover 


a wide area. The liquid crystals were calibrated with 
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respect to the onset of each color. The result is that 
for areas of pulsating or mottled colors, the uncertainty 
of surface temperature and location are much higher then 


forregions in which a distinct color gradient exists. 


Reynolds Number, Re: 
The equation used to calculate Re was: 
Ud 


v 





Re = 


em iicwOunlemuomatrechly proportional to the sqaure 


root of the wall manometer reading, therefore: 





since v iS a property of air, and its uncertainty is 


negligible: 
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Re = 153,000 + 2000 
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miessure COCTTICHENE , Cc: 
The equation used to calculate C was: 








both pressure differences were measured on the same manometer 


bank. For Re = 153,000 and d/r = 0.50: 








oo S 
We ‘ Z , Ze 
p AP, ( 2 | 
== = + 
C AP. AP, 
“(0 
mo S00 |7 mores = OF017 = 0.02 
p p 


C, =-0.99 + 0.02 
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Aree ND A C 
HOT WIRE ANEMOMETER CALIBRATION 


Prior to any experimental runs, the TSI hot wire 
anemometer was calibrated using a controlled air stream and 
a precision manometer. In order to assure accurate calibra- 
tion, the probe was mounted in the wind tunnel traversing 
mechanism and the calibration was done in the wind tunnel 
test section. In this way all leads and connections used 
for experimental runs would be in the circuit for calibration. 

The calibration apparatus provided a controlled air 
flow whose velocity was calculated from the pressure 
differential provided by an alcohol-filled manometer. The 
velocity calculated was plotted against the voltage indicated 
by the hot wire anomometer unit (Figure 20). Calibration 
was done by increasing the calibration apparatus air stream 
velocity so that the alcohol in the manometer would indicate 
even values (in mm). The manometer inclination and readings 
were recorded, along with the voltage indicated by the hot 
wire anemometer. 

The equations used, data taken and calibration curves 
are shown below. Calibration checks were made periodically 
throughout the course of this work, and corrections made to 


account for the aging of the probe. 


Sy) 





EYiommeernoulil) Ss Equation: 


Zia} IM 
U — 
00 p 


where: 
eae otmenmercamevelocity ft/sec 
g. = 32.17 ft/sec’ 
> = air density 1b/ft° 
AP = pressure differential 1b/ft¢ 


The pressure differential from the manometer was calculated 


from the following relations: 








p g 
HA0 
P = ————- _h 
Je 
OS tee, 
and h 0-80 h 
where: 
he = ratio of specific gravity of manometer 
; fluid used, to specific gravity of 
fluid intended for manometer 
C = manometer inclination factor 
h' = height of miniscus in inclined manometer mm 
hence: 
U. = 13.24 vh_ ft/sec 
where: 
h = mm of alcohol 
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Turbulence intensity was determined by: 


c= 4 Ve 
U ye _ ve 
where: 
ae = turbulence intensity 
V = voltage indicated by hot wire anemometer 
3 = rms voltage variation 
Vo = extrapolated value of voltage corresponding 


to zero velocity 
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Voltage sqaured versus the square root of velocity 
for hot wire anemometer calibration. 
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Pie ENIX = D 
WIND TUNNEL VELOCITY PROFILE AND TURBULENCE INTENSITY 


In order to assure that test results would be accurate, 
it was necessary to show that the velocity profile in the 
test section was constant and that turbulence was small 
enough to have no effect on test results. 

The hot wire anemometer was used to probe the test 
section at various locations. The test section coordinate 
axis is shown in Figure 21. The test cylinder was located 
at X = 13.71, Y = -2.33. Velocity profiles were taken at 


the following locations: 


iomemecceomecitemance: A = 17./5 -9.0<Y¥<+9.0 Z = 0.0 
Ova Selo = =2.2<t4 0<2 = 4-0 

PAcnemmcading cage ; eee ee 025-1519), 0 2 = 0.0 

Top of the cylinder: eeeece leo 25<)<79 0229040 


(plate attached) 


icpeaiiGmpomcansor Cyl: X = -]3.7] -10.0<¥<+9.0 Z = 0.0 
(plate at d/r = 5.33) 


All velocity profiles were taken at a tunnel speed of 
98 ft/sec. 

As shown in Figures 22 and 23 the air speed increased 
from 10 to 30 percent on either side of the cylinder in the 
case of the plate at d/r = 5.33, and a like amount over the 
top of the cylinder in the case of the attached plate. The 


velocity at the test section entrance was essentially constant. 
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Turbulence intensity was determined at the test section 


entrance: 
Meare at <>. 0 Z = 0.0 


for three Reynolds numbers; 75,000, 100,000, and 153,000. 
The average value of turbulence intensity for Re = 153,000 
was 0.00239, or about 0.2%. The turbulence for the other 
Reynolds numbers was slightly higher, but of the same order 
of magnitude (see data). The value of turbulence intensity 
was calculated in the same manner as for the hot wire 


anemometer probe calibration described in Appendix C. 
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